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Abstract. Numerical solution of one-dimensional stochastic inte¬ 
gral equations because of the randomness has its own problems, i.e. 
some of them no have analytical solution or finding their analytic 
solution is very difficult. This problem for two-dimensional equations 
is twofold. Thus, finding an efficient way to approximate solutions 
of these equations is an essential requirement. To begin this impor¬ 
tant issue in this paper, we will give an efficient method based on 
Haar wavelet to approximate a solution for the two-dimensional lin¬ 
ear stochastic Volterra integral equation. We also give an example to 
demonstrate the accuracy of the method. 
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1. Introduction 

As we know, two dimensional ordinary integral equations provide an important 
tool for modeling a numerous problems in engineering and science [6, 7]. The second 
kind of two-dimensional integral equations may arise from some problems of nonho- 
mogeneous elasticity and electrostatics [8]. 
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Dobner presented an equivalent formulation of the Dorboux problem as a two- 
dimensional Volterra integral equation [9]. We can also see this kind of equations 
in contact problems for bodies with complex properties [10,11], and in the theory 
of radio wave propagation [12], and in the theory of the elastic problem of axial 
translation of a rigid elliptical disc-inclusion [13], and various physical, mechanical 
and biological problems. Some numerical schemes have been inspected for resolvent 
of two-dimensional ordinary integral equations by several probers. Computational 
complexity of mathematical operations is the most important obstacle for solving 
ordinary integral equations in higher dimensionas. 

The Nystrom method [14], collocation method [15,16,17], Gauss product quadra¬ 
ture rule method [18], Galerkin method [19], using triangular fuctions [20, 21], Legen- 
der polynomial method [22], differential transform method [23], meshless method [24], 
Bernstein polynomials method [25] and Haar wavelet method [26]. This paper is hrst 
focused on proposing a generic framework for numerical solution of two-dimensional 
ordinary linear Volterra integral equations of second kind. The use of the Haar wavelet 
for the numerical solution of linear integral equations has previously been discussed 
in [1] and references therein. The paper [1] should be considered as a logical contin¬ 
uation of the papers [2 — 4]. In [2] a new numerical method based on Haar wavelet 
is introduced for solution of nonlinear one-dimensional Fredholm and Volterra inte¬ 
gral equations. In [3] the Haar wavelet method [2] is extended to numerical solution 
of integro-differential equation. In [4] the Haar wavelet method [2, 3] is improved in 
terms of efficiency by introducing one-dimentional Haar wavelet approximation of the 
kernel function. The method [1] is fundamentally different from the other numerical 
methods based on Haar wavelet for the numerical solution of integral equations as it 
approximates kernel function using Haar wavelet. 

The general hyperbolic differential equation is dehned as [9] 

(1.1) u^y = r{x,y,u,Ux,Uy) , {x,y) e B 

u^^ = h{x,y) , BCdB, 

where the domain B and the subset B of the border dB are chosen according to the 
different initial value problems. It’s easy to show that the integral form of (1.1) is 
given by two-dimentional Volterra integral equation 

n x 

Ki{x,y,s,t)g{s,t)dsdt. 
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Similarly, if we import statistical noise in to (1.1), we can obtain two-dimensional 
linear stochastic Volterra integral equation of the second kind, i.e. 


( 1 . 2 ) 


9{.x,y) = f{x,y) + 

+ 





'0 JO 
ry 



0 Jo 


Ki{x, y, s, t)g{s, t)dsdt 

K2 {x, y, s, t)g{s, t)dB{s)dB{t) 

(x, y) e [0,1] X [0,1] , s^x <t^ y. 


where the kernels Ki{x,y, s,t) and K2{x,y, s,t) in (1,2) are known functions and 
f{x, y) is also a known function whereas g{x, y) is unknown function and is called the 
solution of two-dimensional stochastic integral equation. The condition s ^ x < t ^ y 
is necessary for adaptability to the hltration {Fy, 0 < t < 1} where Ft = a{B{s)-, 0 < 
s < 1}. 


Lemma l.Put s) = K{x, y, s, t)g{s, t). Let 0 be a function in T^([0,1]^). Then 
there exists a sequence 0^ of off-diagonal step functions such that 


lim 

n^oo 



0(f,s) - 0n(t,s) 


dtds 


0 . 


Definition 1. Let 0 G T^([0,1]^). Then the double Wiener-Ifb integral of 0 is 
dehned as 


h pb 



r*b pb 


s)dB{t)dB{s) = lim 


(pnit, s)dB{t)dB{s) in L^(f2). 


2. Haar Wavelets 


A wavelet family (0j,j {y))j£N,i£Z is an orthonormal subfamily of the Hilbert space 
L‘^{R) with the property that all function in the wavelet family are generated from a 
hxed function 0 called mother wavelet through dilations and translations.The wavelet 
family satishes the following relation 


^jn {y) = ■ 


For Haar wavelet family on the interval [0,1) we have: 


hi{y) 


1, forje [0,1) 
0, otherwise, 
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hi{y) = 


1 

- 1 , 

0 , 


forye [a, 13) 
forye [/3,y) 

otherwise, i = 2,3,,..., 


a = 


n 

m’ 


m = 




(n + 0.5) 
m 


7 = 


(n + 1) _ 
m ’ 


£ = 0,1,..., n = 0,1,..., m — 1. 


The integer £ indicats the level of the wavelet and n is the translation parameter. 
Any square integrable function fiy) dehned on [0,1) can be expressed as follows: 

OO 

i=l 

where a* are real constants. 

For approximation aim we consider a maximum value L of the integer £, level of the 
Haar wavelet in the above dehnition. The integer L is then called maximum level of 
resolution. We also dehne integer M = 2^. Hence for any square integrable function 
f{y) we have a hnite sum of Haar wavelets as follows: 

2M 

fiy) ~ ^aiHy)- 

i=l 

The following notation is introduced [1]: 

ry 

(2.1) PiAy) = / hi{u)du, 

Jo 

where by the dehnition of Haar wavelet equation (2.1) reduce to 


Pi Ay) 


y-a, for y e [a, (3) 

1-y, fory^[l3,y) 

0, elsewhere. 


We can also the following stochastic notation introduce 

(2.2) hiAy) = [ hi{u)dB{u), 

Jo 

where equation (2.2) can be evaluated similarly by the dehnition of Haar wavelet and 
is given as follow: 
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r B{y)-B{a), 

%i{y) = \ 2B{li)-B{a)-B{y), 

I 0 ^ 


for y e [a, 13) 
for y e [13, 7 ) 
elsewhere. 


3. Numerical method 

In this section, proposed nnmerical method [1] will be discussed for two-dimensional 
linear stochastic Volterra integral equation of the second kind. In the hrst subsec¬ 
tion, we state some results for efficient evalution of two-dimensional Haar wavelet 
approximations. In the second subsection, we apply these results for hnding numeri¬ 
cal solutions equation ( 1 . 1 ). 

For Haar wavelet approximation of a function f{x,y) of two real variables x and y, 
we assume that the domain 0 < x, y < 1 is divided into a grid of size 2M x 2N using 
the following collocation points 


(3,1) 

Xm= 2 ^ ,m = 1 , 2 , ..., 2 M, 

(3.2) 

yn= ,n = l,2,...,2N. 


3.1 Two-dimensional Haar wavelet system 

A real-valued function G{x, y) of two real variables x and y can be approximated 
using two-dimensional Haar wavelets basis as [1,27]: 

2M 2N 

(3.3) G{x,y)^ EE ^p,qhp{x)hq{y) . 

p=l q=l 

In order to calculate the unknown coefficients the collocation points dehned 

in Eqs. (3.1) and (3.2) are substituted in Eq. (3.3).Hence, we obtain the following 
2M X 2N linear system with unknowns 

2M 2N 

(3.4) G{xm,yn) = 'Y^'Y^hp^qhp{xra)hq{yn),m = 1,2,...,2M, n = l,2,...,2N. 

p=l q=l 

The solution of system (3.4) can be calculated from the following theorem. 
Theorem 2. The solution of the system (3.4) is given below: 
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2M 2N 


bi,i = 


2M X 2N 


EE yn)i 


p=l q=l 


/3i 2N 


71 2N 


bi,i = 


bi,j = 


bi,j — 


Pi X 2N 
1 

2M X p2 


^^G{xm,yn)- ^ '^G{x„,,yn)\ , i = 2,3,...,2M, 

\p=ai q=l p=/3i + l q=l 

' 2M P 2 2M 72 

'^^G{Xm,yn) ^ G{Xm,yn)\ , j = 2, 3, ... ,2N, 

iP=l 9=02 P =1 9=/32+1 


Pi X P2 \ “ 


( /3i P 2 Pi 72 71 P 2 

EE Vn) E E Vn) E Ec(-™.!/. 

p= 0 i 9=02 P=oi q=p 2 +l p=/3l+l 9=“2 


71 72 


- ^ ^ G'(a;m,|/„) j , i = 2,3, ...,2M , j = 2, 3,..., 2iV, 

p=/3i+l 9=/32+1 


where 


(3.5) 


and similarly, 


(3.6) 


Proof. See [2], 


cii — Pi(ii'i ~ 1) + 1) 
A = Pi(o'i — 1) + Y, 
7i = Pio-i, 

(Tl = i - Ti, 

= 2 l-'°§ 2 (*-i)J 


0(2 = P2{(T2 - 1 ) + 1 , 
^2 = P 2 ( 0'2 — 1 ) + Y, 
72 = P2<72, 

n - 21 V 

P 2 j-2 ) 

<72 = j - 72, 
r 2 = 2 L'°g 20 '-i)J 


□ 
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Consider a function G{x, y, s, t) of four variables x, y, s and t. Suppose G{x, y, s, t) 
is approximated using two-dimensional Haar wavelet as follows [1]: 

2M 2N 

(3.7) G{x,y,s,t)^ EE bp^q{x, y')hp[s')hq(t'). 

p=l q=l 


Substituting the collocation points 

i-0.5 


Si — 


and 


2M 
j - 0.5 

ti = - - 

^ 2N 


, j = l,2,...,2iV, 


we obtain the linear system 
(3.8) 

2M 2N 

G{x,y,Si,tj) ^ EE bp,q {x,y)hp{si)hq{tj) , i = l,2,...,2M , j = 1, 2,2iV. 

p=l q=l 


Corollary 1. The solution of the system (3.8) for any value of x.y G [0,1] is given 
as follows [ 1 ]: 


2M 2N 


bi,i{x,y) = 


Pi 2N 


2M X 2N 


^^G{x,y,Sp,tq), 


p=l 9=1 

71 2N 


bi,i{x,y) = 
bi,j{x,y) = 
bi,j{x,y) = 


.,x 2 iVlEE G{x,y,Sp,tq) - E E G{x,y,Sp,tq) j ,7 = 2,3, ...,2M, 

Vp=oi 9=1 p=/3i+l 9=1 

^ 2M 02 2M 72 

^ ^ ^ ^ G{x, y, Sp, tq) 'y ^ y ^ G{x, y, Sp, tq) | , j = 2, 3,..., 2iV, 

vP=l <?=02 p=l 9=/?2 + l 


2M X P‘2 
1 


Pi X p2 


/ /3i ^2 /3i 72 

EE G{x,y,Sp,tq) - 

\p=Qi9=^2 P=ai 9=P2+1 


71 72 


“qj ) 


7i P2 

E E G{x,y,Sp,tq) + 5 : 5 : G(x,v,s„t, 

p=l3\ + lq=a 2 P=/3i-l-l 9=/32-I-1 

7 = 2,3,...,2M , j=2,3,...,2iV, 


where oi, /3i, 71 and pi are dehned as in Eq. (3.5) and 0 : 2 , /d 2,72 and p 2 are dehned as 
in Eq. (3.6). 

Corollary 2. Suppose a function G{x,y) of two variables x and y is approximated 
using Haar wavelet approximation given in Eq. (3.3). Suppose further that G{x,y) 
is known at collocation points {xm,ym), tu = 1,2, ...,2M,n = l,2,...,2iV. Then the 
approximate value of the function G{x, y) at any other point of the domain can be 



calculated as follows [1]: 


2M 2N 


G{x,y) = 


G{xm,ym)hi{x)hi{y) 


2M ^ / /3i 27V 


E 


2M X 2N 

p=l q=l 

71 27V \ 

-.,x2iVlEE G(Xm) ym) E E G{xm,ym) I hi{x)hi{y) 

7=1 \p=ai q=l p=/3l+l q=l / 

27V ^ / 27U /32 27U 72 \ 

EE G 1/m) E E G(x^,ym) j /ii(a;)/ij(i/) 

j = l Vp=l 9=^2 «?=/32 + 1 / 

2M 2Ar / ^2 /3i 72 

+ ym) - F{xm,ym) 

i=l j=l ^ ^ \p=oi q=a2 p=ai g=/32+l 

71 /32 71 72 \ 

- E E Gi^X'fji^ yrfi) H“ E E G{xm,ym) j hi{x)hj{y), 

p=/3i+lg=02 p=/3i+l (j=/32+1 / 


where ai, /3i, 71 and pi are defined as in Eq. (3.5) and 0:2, P 2,12 and p 2 are defined as 
in Eq. (3.6). 

3.2 Two-dimensional linear stochastic Volterra integral equation 


Consider the two-dimensional linear stochastic Volterra integral equation (1.2). As¬ 
sume that the function K{x, y, s, t)g{s, t) is approximated using two-dimensional Haar 
wavelet as follows: 

2M 2N 

(3.9) Ki{x,y,s,t)g{s,t) ^ EE hi^j{x,y)hi{s)hj{t). 

7=1 j=l 


2M 2N 

( 3 . 10 ) K 2 {x,y,s,t)g{s,t) ~ EE Ci^j{x,y)hi{s)hj{t). 

7=1 j=i 

With this approximation Eq. (1.2) can be writen as follows: 

n x 2717 27V 

Y] Y y)hi{s)hj{t)dsdt 
7=1 j=i 



2M 27V 

EE Cij {x,y)hi{s)hj (t)dB {s)dB (t). 

7=1 j=l 


Eq. (3.11) can be written in a more compact form using the notations introduced in 
equations (2.1) and (2.2) and is given as follows: 

2M 27V 2M 2N 

9{x,y) = f{x,y) + EE bi,j{x,y)pi^i{x)pj^i{y) + EE Ci,j{x,y)qi^i{x)qj,i{y). 

2=1 j = l 2 = 1 j=l 
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Substituting the collocation points given in (3.1) and (3.2), we obtain the following 
system of equations: 

2M 27V 

di^m^y-n) fi^m^yn) 3“ ^ ^ ^ l/n)P7,l (^m)Pj,l (?/n) 

7=1 j=l 


2M 2N 

Ci,j{x^,yn)qi,l{Xm)qj,l{yn)- 

7=1 j=l 

Now , i = l,2, ...,2M , j = 1, 2,..., 2iV and similarly C 7 j , i = 1, 2,..., 2M , j = 
1,2, ...,2iV can be replaced with their expressions given in Corollary 1 and the fol¬ 
lowing system of equations is obtained: 

, . , . 2M 2N 

(3.12) gix^,yn) = f{Xm,yn) + l/n, tq)g{Sp, tq) + 

p=l q=l 


2M , . , . / pi 2N 

g7,l(a^m)Pl,l(|/n) 


7=2 


Pi X 2N 


2N 

E 

i=2 


Pl,l{Xm)Pj,l{yn) 

2M X P2 


2M 2N 

EE 


Pi,l{Xm)Pj,l{yn) 
Pi X p2 


^ ^ ^ ^ ^1 {Xmi Pni Sp, tq^gi^Sp, tq) 

p=a\ q=l 

71 27V 

^ ^ ^ ^ KliXrni ym Spi tq) g{Sp,tq) | -|- 

P=/3i + l 9=1 

2M P 2 

^ ^ ^ ^ ^l{Xmiyni Sp,tq) g{Sp,tq) 
p=lp=a2 

2M 72 

'y ^ ^ ^ Kl{Xm)yn) Spitq)g(^Sp^tq) | -|- 

p=l q=P2 + i 

01 02 

y ^ y ^ ^liXmiyni Spitq)g{^Sp,tq) 

p=oi (7=02 

7l 02 


i—2 j—2 
Pi 72 

y ^ ^ ^ Kii^Xmtyni Sp,tq) g{Sp,tq) ^ ^ Kl{Xmi yni Sp, tq)g{Sp, tq)-\- 

p=ai q =02 + l p= 0 i + lq=a 2 

71 72 \ 

y ^ ^ ^ Kl{Xm,yn, Sp,tq')g(^Sp,tq) j -|- 

p=pl + l q=p2 + l / 

. ^ . 2M 2N 

qi,i[x„i)qip[yn) 


2M X 2N 


EE ^2iXmi yni ^pi tq)gi,^pi 7)”^ 


p=l (j=l 


2M , . , . / 0x 2N 

yy g7,l(a7m)gi,l(l/n) 


7=2 


Pi X 2iV 


EE K^2 [Xmi yni Spi tq)g(ySp^ tq^ 

p=ai q=l 

7i 27V 

E E ^2iXmiyni Spitq) g(^Sp,tq 

p=0i + l q=l 
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2N , . , ^ /2M h 


3=2. 


2M X p2 


^2iXmi Vni ^pitq) 

p=\ p=a2 


2M 72 

E E K^2{Xm)yn) Sp)iq) 9{,^p)^q) | E 

p=l q=/32 + l 


2M 2N f ^ f ■,/ Pi h 

qi,i\Xm)qj,i[u 


pi X p2 


EE ^2iXmiyni Spitq) y{Sp,tq)- 


\p=ai q=a2 


i=2 j=2 

Pi 72 71 p2 

E E K^2{Xmiyni Spitq) di^Sp^tq) E E Vni ^pitq) 9{,^pi^q)~^ 


p=ai g=^2 + l 


p=/5i + l q=<^2 
7l 72 


E E ^2iXmiyni Spitq) g(^Sp,t 

p=Pl + l (J=/32 + 1 


q/ ’ 


m = l,2,...,2M , n = 1, 2,..., 2iV. 


Eq. (3.12) represents 2M x 2iV system which can be solved using either prevalent 
methods for solving linear systems. The solution of this system gives values of g{x, y) 
at the collocation points. The values of g{x, y) at points other than collocation points 
can be calculated using Corollary 2. 


4. Numerical Example 


In this section, the numerical example is given to demonstrate the applicability and 
accuracy of our method. Consider the following linear 2D stochastic Volterra integral 
equation of second kind: 

n x ny PX 

{x+y+t—s)u{s,t)dsdt+ / {x+y+t+s)u{s,t)dB{s)dB{t) 

Jo Jo 


where 

f{x, y) = x + y- —xy{x^ + Ax‘^y + Axy"^ + y^). 

The solutions mean together conhdence interval at the collocation points for the 
present method for 1000 iterative of system (3.12) is shown in Table 1. In Figs. 1 — 4, 
three-dimensional graphs of the approximate solution for various values of level L are 
shown. 


Table l:The solutions mean together conhdence interval for above example 
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J 

M 

2 M 


u{x,y) 

L 

Confidence 

Interval U 

0 

1 

2 

(0.25,0.75) 

1.9951 

1.9951 

1.9951 




(0.125,0.375) 

1.06717 

1.06715 

1.06719 

1 

2 

4 

(0.375,0.875) 

2.43481 

2.43194 

2.43769 




(0.625,0.875) 

2.72591 

2.72036 

2.73147 




(0.0625,0.4375) 

1.0498 

1.0498 

1.04981 

2 

4 

8 

(0.3125,0.6875) 

2.15292 

2.15266 

2.15318 




(0.8125,0.9375) 

2.68057 

2.67707 

2.68407 




(0.03125,0.71875) 

1.51492 

1.51488 

1.51495 

3 

8 

16 

(0.40625,0.53125) 

2.25846 

2.25532 

2.26159 




(0.78125,0.96875) 

2.67147 

2.67139 

2.67156 




(0.015625,0.609375) 

1.26132 

1.26129 

1.26135 

4 

16 

32 

(0.296875,0.796875) 

2.25999 

2.25891 

2.26107 




(0.859375,0.984375) 

2.58876 

2.58870 

2.58881 




4 


Figure 1. Plot of approximate solution of level L = 1 for test example. 



Figure 2. Plot of approximate solution of level L = 2 for test example. 
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Figure 3. Plot of approximate solution of level L = 3 for test example. 



Figure 4. Plot of approximate solution of level L = 4 for test example. 
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5. Conclusion 

As mentioned above, numerical solution of two-dimensional stochastic integral 
equations because of the randomness is very difficult or sometimes impossible. In 
this paper, we have successfully developed Haar wavelets numerical method for ap¬ 
proximate a solution of two-dimensional linear stochastic Volterra integral equations. 
The example conhrm that the method is considerably fast and highly accurate as 
sometimes lead to exact solution. Although, theoretically for getting higher accuracy 
we can set the method with larger values of M and N and also larger of the degree 
of approximation, p and q, but it leads to solving MN linear systems of size pq x pq, 
that have its difficulties. The method can be improved to be more accurate by using 
other numerical methods. Mathematica has been used for computations in this paper. 
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